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Abstract. It is shown that the second term in the asymptotic expan- 
sion as f ^ of the trace of the semigroup of symmetric stable processes 
(fractional powers of the Laplacian) of order a, for any < a < 2, in 
Lipschitz domains is given by the surface area of the boundary of the 
domain. This brings the asymptotics for the trace of stable processes in 
domains of Euclidean space on par with those of Brownian motion (the 
Laplacian), as far as boundary smoothness is concerned. 

1. Introduction and statement of main result 

Let Xt be a symmetric a-stable process in W^, a G (0,2]. This is a pro- 
cess with independent and stationary increments and characteristic function 
^Ogi^Xt ^ e-t|€r , ^ G M^, t > 0. By p(t, x, y) = pt{x - y) we will denote the 
transition density of this process starting at the point x. That is, 



p-{XteB)= [ p{t,x,y)dy. 
Jb 



Since the transition density is obtained from the characteristic function by 
the inverse Fourier transform, it follows trivially that pt{x) is a radial sym- 
metric decreasing function and that 

pj(x) = r'^/>i(t"^/"x) < t-'^/>i(o), t>o, xeR'^. (1.1) 

We also have (see (1.2) [l]) 

where uJa = r{d/2) surface area of the unit sphere in W^. 

Let D C be an open nonempty set and denote by Tn = 'mf{t > : Xt ^ 
D} the first exit time of Xt from D. The transition density poit, x, y) of the 
process killed while exiting a domain D (a-stable heat kernel) is defined by 

P'^iXt eA,TD>t)= [ pD{t,x,y)dy. (1.2) 
J A 

This subprobabilistic density satisfies 

PD{t, X, y) = p{t, X, y) - rnit, x, y), (1.3) 
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where 

rD{t,x,y) = E^TD <t;p{t-TD,X{TD),y)). (1.4) 

By {Pt^}t>o we denote the semigroup on L?'{D) of Xt killed upon exiting 
D. That is, for any t > and / G L^{D) we define 

Pt''f{x) = E^TD>t;f{Xt))= [ pD{t,x,y)f{y)dy, x G D. 

Jd 

The study of the spectral properties of the a-stable heat semigroup {PP}t>o 
has been the subject of many papers in recent years see e.g. [1], [2], [9], [iU], 
[12]. Whenever D is bounded (or of finite volume), the operator Pf^ maps 
L'^{D) into L°°{D) for every t > 0. This follows from (1.1), (1.4), and the 
general theory of heat semigroups as described in [11]. In fact, it follows 
from [11] that there exists an orthonormal basis of eigenfunctions {(pn}^=i 
for L'^{D) and corresponding eigenvalues {AnjJ^iof the generator of the 
semigroup {PP}t>o satisfying 

< Ai < A2 < As < . . . , 

with An — > 00 as n ^ 00. That is, the pair {(pm^n} satisfies 

PrMx) = e-^"Vn(x), x£D, t>0. 

Under such assumptions we have 

00 

poit, x,y) = Y, e-^"^^n{x)ipn{y)- (1.5) 

n=l 

The trace of the a-stable heat kernel on D (often referred to as the par- 
tition function of D) is defined by 

ZD{t) = I PD{t,x,x)dx. (1.6) 
Jd 

Because of (1.5), we can rewrite (1.6) as 

00 „ 00 

Zoit) = Y^e-'-' / ^lix)dx = J^e-^"*. (1.7) 

n=i 

It is shown in [ !] that for any open set D C M*^ of finite volume 

\iint'^/'^ZD{t) = Ci\D\, (1.8) 

where Ci = uJd^{d/a){2'K)~'^a~^ . This result is proved in [l] under the 
assumption that dD has zero Lebesgue measure. As observed in Remark 
2.2 of [ -], the result in fact holds for all open sets of finite volume. 

As is well known, the asymptotic behavior of the partition function as 
t — >■ implies Weyl's formula on the growth of the number of eigenvalues. 
Indeed, if we let A^(A) be the number of eigenvalues {Aj} which do not exceed 
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A, it follows from (1.8) and the classical Karamata tauberian theorem (see 
[2] for full details) that 

^(^)~ WTT^^'^"' asA-oo. (1.9) 
1 [a/ a + 1) 

This is the analogue for stable processes of the celebrated Weyl's asymptotic 
formula for the eigenvalues of the Laplacian. 

The asymptotics for the trace of the heat kernel when a = 2 (the case of 
the Laplacian with Dirichlet boundary condition in a domain of R"^), have 
been extensively studied by many authors. In particular, van den Berg [ ] 
proved that under an assumption of i?-smoothness of the boundary (that 
is dD satisfies uniform outer and inner ball condition with radius ii), when 
a = 2, 



< _£J ! f -> 



(1.10) 

When the domain has boundaries the result 



Znit) = (47rt)-^/2 J^i^i _ v^^d-i(g^) ^ ^(^1/2^^ , t ^ 0, (1.11) 

was proved by Brossard and Carmona in [(>]. R. Brown subsequently ex- 
tended (1.11) to Lipschitz domains in ["]. We refer the reader to [ ] and [7] 
for more on the literature and history of these type of asymptotic results. 

In [2], the exact analogue of (1.10) was proved for any stable processes 
of order < a < 2. Our goal in this paper is to obtain the second term in 
the asymptotics of Zj:){t) under assumption that D is a bounded Lipschitz 
domain, in complete analogy to the R. Brown result, [7]. 

Theorem 1.1. Let D C M^, d > 2, be a bounded Lipschitz domain. Let 
\D\ denote d- dimensional Lehesgue measure of D and 7i'^~^{dD) denote the 
{d — 1)- dimensional HausdorjJ measure of dD. For any < a < 2, the 
partition function of the symmetric a-stable process in D satisfies 

f^/'^Zoit) = Ci\D\ - C2n'^-^{dD)t^''' + o(t^/°), (1.12) 



where 



and 



Here, 



ujdT{d/a) 

Ci=Pi(0) = ^^^ (1.13) 

POO 

C2= / rH(l,(xi,0,...,0),(xi,0,...,0))fixi. (1.14) 

^0 



H = {{xi,...,Xd) G M"* : xi > 0} = 
is the upper half-space o/M'^ and 

rH{t,x,y) = E'^{th <t;p{t-TH,X{TH),y)), (1.15) 
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is as in (1.4)- 

It is now well known (see [15], [17], [18], [19], [20] and references therein) 
that in the case of the Brownian motion (the Laplacian), the second term 
in Weyl's asymptotics for A^(A) is also given by the surface area of the 
domain-at least in the case of smooth domains. In fact, 

N{\) = Ci\D\ X^''^ - C2\dD\\^ + o{\^), (1.16) 

where Ci and C2 are explicit constants depending only on d. 

The arguments used to obtain such results employ tools from the theory 
of the wave equation. An interesting and challenging problem is to develop 
the wave techniques in the case of the fractional Laplacian to obtain simi- 
lar results for the counting function of stable processes. To the best of our 
knowledge the "wave group" corresponding to these operators has not been 
studied before. An alternative approach would be to find a probabilistic 
(heat equation) proof for (1.16) and then try to adapt such arguments to 
stable processes. A success with either approach is likely to lead to applica- 
tions that will be of independent interest. 

2. Preliminaries 

In this section we present several facts concerning symmetric a-stable 
processes and recall several geometric properties of Lipschitz domains which 
will be needed in the proof of our main result. Theorem 1.1. These geometric 
facts, and notation, for Lipschitz domains are standard and follow [7]. 

The ball in with center at x and radius r, {y S M'^ : |x— y| < r}, will be 
denoted by B{x^ r). We will use 5d{x) to denote the distance from the point 
X to the boundary, dD^ of D. That is, 5d{x) = dist(x,OD). Throughout 
the paper, we will use c to denote positive constants that depend (unless 
otherwise explicitly stated) only on d and a but whose value may change 
from line to line. 

The Levy measure of the stable processes Xt will be denoted by v. Its 
density, which we will just write as J^(x), is given by 

= j^, (2.1) 

where Ad,'^ = T{{d - 7)/2)/(2T7r'^/2|r(7/2)|). We will need the following 
bound on the transition probabilities of the process Xt which can be found 
in [23]. For ah x,y and t > 0, 

The scaling properties from (1.1) of pt{x) are inherited by the kernels po 
and r£). Namely, 



1 / X y 

PDit,x,y) = jjj^PD/t^/" (1.^ 



/a ' J ' 
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and 

rB(t,x,y) = ^r^/,i/. (l,^,^). (2.3) 

Also, both pd and td are symmetric. That is, pD{t,x,y) = pD{t,y,x) and 
r£){t,x,y) = r£,{t,y,x). The Green function for the process in the open 
set L) C will be denoted by G£){x,y). In fact, this can be written in 
terms of the transition probabilities as 



GD{x,y) 



/•oo 

/ PDit,x,y)dt, x,ye 
Jo 



We have the following estimates for r^. 
Lemma 2.1. Let D CW^ be an open set. For any x,y €z D we have 

'■"''•-•'''^'^(^'^ <^-^' 

This Lemma follows from [ ], Lemma 2.1. 

Lemma 2.2. Let D C M.'^ be an open nonempty set. Fix e > 0. For any 

y £ D, and x,z £ D such that 6d{x) > e, 6d{z) > e, we have 

|rD(l,x,y) - rD{l,z,y)\ < c{e)\x - z\. 
Here c(e) depends on e, d, a. 
Proof. Recall that 

where gt : (0, oo) (0, oo) is the density of the a-stable subordinator whose 
Laplace transform is given by e~^^gt{u) du = e^*'^"''^. 

For the rest of the proof of this Lemma, let us denote pi by p^^ in order 
to stress the dependence on the dimension d. Differentiating in (2.5), we get 
for any x G M*^, 

l^(^) = -2^1 / (4.,)I.+2)/2 (^) 51 (^) = -2x,^ir'\x), 

where x = (x,0,0) G 

Since for any dimension d we have (see (2.2)), p^f^ {x) < c|x|~'^~", we get 



^''^ (X) 



dxi 



< c\xi\\x\-'^-^-'' < clxl-'^-^-". 
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Set ei = (1,0, ...,0) G M"*. Since the function pf^ is radial, p'i\x) 



p\ (|2;|ei). Hence, the mean-value theorem gives 



Ip'iHx) - pf {z)\ = Ip'f' {\x\ei) - p'f' {\z\ei 



id), 



dp\ 



(d) 



dxi 

< c||j;| A \z\\~ 



- \z\ei\ 



\x\ — \z\ 

d-l-a\ 



\X\ — \Z\ 



where ^ = (Ci,0, ... ,0) G M'^ is the point between \x\ei and \z\ei. By the 
scaling of pi , 

rr,(l,x,y) = Ey {td < 1; p{l - td , X (td) , x)) 



Ey[TD< 1; 



1 



[d) ( X- X{td) \ 



It follows that 

\rD{l,x,y) - rz)(l,z,y)| 
1 



= Ey iyTD < 1; 
< cEy I TD < 1; 



(1 - tdY/' 

1 



id) 



nrp) \ _ ^id) ( z-X{td) \ 
(l-rz5)V-; 



X-XiTD)\ ^^ \z-XiTD)\ \ 

(1 - td)'^/° V (1 - td)^/" (l-T£))l/"y 
x-X{td)\ \z-X{td)\ 



-d-l-a 



(1 - rz))i/" (1 - rfl)V" 

< cEy{TD < 1; (1 - TD))(fe(x) A 

< c{e)\x — z\, 



X — z 



where we used our assumption that both doix) and 5d{z) are larger than 
e. □ 

As an immediate corollary of this lemma we obtain the following diagonal 
estimate for r^. 

Lemma 2.3. Let D C M.'^ be an open nonempty set. Fix e > 0. For any 
x,z £ D such that 5d{x) > e, Sd{z) > e, we have 

\rD{l,x,x) - rD{l,z,z)\ < c{e)\x - z\. 

Here c{e) depends on e, d, a. 

Proof. By the fact that r/)(l,x,2:) = r£){\,z^x) and Lemma 2.2 we get 

\rD{l,x,x) - rD{l,z,z)\ < |rD(l, x, x) - rD(l, 2;, x)| + |rz)(l, x, z) - r£)(l, 

< c{e)\x — z\. 

□ 
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In the sequel we will need several facts concerning the a-harmonic mea- 
sure P''^{X[td) G •) of an open nonempty set D <Z W^. We say that an 
open set D <zW^ satisfies the outer cone condition if there exist constants 
T] = rj{D), Ro = Ro{D) and a cone C = {x = {xi,...,Xd) e M.'^ : < 
Xd, Wixi, . . . ,Xd-i)\\ < rjXd} such that for every Q G dD, there is a cone Cq 
with vertex Q, isometric with C and satisfying Cq n B{Q,Rq) C V^. It is 
well known that if D is an open set satisfying the outer cone condition then 

P'-'iTD < oo, X{td) edD) = 0, xe D. 

This fact is proved in [ ], Lemma 6, for bounded open sets with the outer 
cone condition. However, by the same arguments as in the proof of Lemma 
2.10 in [1(3], it holds also for unbounded open sets. 

In the sequel we will also use the Ikeda-Watanabe formula for the space- 
time a-harmonic measure in terms of the Green function GD{x-,y) (or the 
transition density p£i{t,x,y)) and the Levy measure of the process. This 
formula is from [ !2], Theorem 2.4. 

Proposition 2.4. Let D be an open nonempty set and A a Borel set such 
that A C D'^ \ dD. Assume that < ti < t2 < oo, x £ D. Then we have 

P'-'iti <TD <t2, X{td) £ A) = f [ pD{s,x,y)ds [ v{y-z)dzdy. 

Jd Jti J A 

(2.6) 

By letting ^2 — > oo in (2.6) we see that this formula also holds for t2 = oo. 
With this and the fact that 

/•oo 

GD{x,y) = / pD{t,x,y)dt 
Jo 

we have 

Corollary 2.5. Let D be an open nonempty set and A a Borel set such that 
A C D'^ \ dD, X € D. Then we have 

P^'iTD <^,X{td) £ A) = / GDix,y) / u{y-z)dzdy. 

Jd J a 

We also recall here that 

Kd{x,z)= [ GD{x,y)u{y-z)dy (2.7) 
Jd 

denotes the Poisson kernel for the a-stable symmetric process and the open 
set D. 

As in the proof in [7], we need to divide the domain D into a good and 
a bad set. We recall several geometric facts about Lipschitz domains, most 
which come from [ ]. 

Definition 2.6. Let e, r > 0. We say that G C dD is {£,r)-good if for each 
point p £ G, the unit inner normal ^{p) exists and 

B{p,r)ndD C {x : \{x - p) ■ u{p)\ < e\x - p\} . (2.8) 
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Using this definition we can construct a good subset of the points near 
the boundary. 

g=[Jrrip,e), (2.9) 
peG 

where Tr{p,e) is a truncated cone 

Trip, = : (x - p) ■ v{p) > \/l - e^\x - p| | n B{p, r) (2.10) 

Figure 1 shows one of the cones r,.(j), e) together with the boundary of D 
contained m. {x : \ {x — p) ■ i'{p)\ < e\x — p\}. 




Figure 1 . Cones at a good point p. 



Lemma 2.7. Let < e < 1/2, r > and suppose that G is a measurable 
{e,r)-good subset ofdD, the boundary of a Lipschitz domain. There exists 
S(){dD, G) such that for all s < sq 

\{xe D:6d{x) <s}\g\< siW^-^dD \ G) + e(3 + n'^-\dD))]. (2.11) 

This lemma states that the measure of the set of the bad points near 
the boundary is small. The proof is essentially the same as the proof as 
Proposition 1.3 in [ . ]. We present the proof here due to the fact the set Q 
is constructed from narrower cones than those in [7]. 

Proof. It is enough (see proof of [ , Proposition 1.3]) to show that 

\{x£D: 5d{x) <s}ng\> siT-C^-^G) - 6(2 +n'^-\G))). (2.12) 
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Choose i>i,i'2, . . . ,1'N £ S'^^^ and closed disjoint sets Fi, F2, . . . , F/v C G 
satisfying H'^^^{G \ IJili ^i) < ^1 diam(Fj) < r and \v{p) — Ui\ < e, for any 

Note that p + pvi G Tr{p.,e) for any p < r and p & Fi. Indeed e < 
1(0, 1) — (e, \/l — e'^)\, where the last distance is the distance between iy{p) 
and the arbitrary unit vector ppi, where pi is on the boundary of Too(p, e). 

We claim that p ^ p + pvi is injective when p ^ Fi and < p < r. 
Suppose that p — q = {p — for p,q £ Fi and p,(T £ (0, r). Since p £ G 
and ll? — (?| < r 

\p - a\ = \ {p - a)ui ■ i^il = \{p - q) ■ Ui\ 

<\[p-q)- u{p)\ + \{p-q)-{iyi- iy{p))\ (2.13) 
< e\p — q\ + e\p — q\ = 2e\p — a\. 

This is a contradiction for e < 1/2. 

Now we can estimate the Lebesgue measure of the set {x £ D : 5d{x) < s}n 
Q. Let 7 denote the minimum of 2 dist(i<i, -Fj) for i ^ j and sq = min(r, 7). 
For any s < sq 

N 

\{x£ D : 5d{x) <s}ng\>^\{p + tiyi:p£Fi,0<t< s}\ 

(2-14) 

> s(l-e)^?^'^-i(Fi) > sin^-\G)-e{2 + n'^-\G))), 
1=1 

where the second inequality follows from the co-area formula [L-j, Theorem 
3.2.3] applied to the map {p, t) ^ p + tvi. This proves (2.12). □ 

The existence of a set G is established in [7, Section 4]. We recall the 
following from [ ], Section 4, page 897. 

Lemma 2.8. For arbitrary e > there exists r such that an {e,r)-good set 
G exists and 

n'^-\dD\G) <e. (2.15) 
This and Lemma 2.7 imply that 

\{x£D: 5d{x) <s}\g\< se{A + n'^-\dD)). (2.16) 

For arbitrary < e < 1/4, let G be the (e,r)-good set from the above 
Lemma. We construct a good set Q using this particular set G. For any 
point x\nQ there is a point on the boundary p{x) such that x £ Tr{p{x)^e). 

We define inner and outer cones as follows 

Ir{p{x)) = {y:{y- p{x)) ■ u{p{x)) > e\y - p{x)\} H B{p{x),r), (2.17) 

Uripix)) = {y:{y- p{x)) ■ u{p{x)) < -e\y - p{x)\} n B{p{x), r). (2.18) 
For e < 1/4 and x G ^ we have 

Tr{p{x),e) C Ir{p{x)) CDC U^{p{x)). (2.19) 
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Now our aim is to show that there exists a half-space H*{x) such that 
xeH*{x), 6H*i,){x) = 6D{x), Ir{p{x))cH*{x)cU^{p{x)). (2.20) 




Figure 2. A halfspace preserving the distance from x to the 
boundary of D. 

Let H^{x) be the halfspace containing the inner cone Ir{p{x) with the 
boundary dH^{x) containing p{x) and perpendicular to x — p{x). Con- 
sider an arbitrary 2-dimensional plane containing the axis of the inner cone 

Ir{p{x)). (The projections of Ur{p{x)), Ir{p{x)), Tr{p{x),e) onto this plane 

> 

are shown on Figure 1.) Let p{x)pi be a unit vector lying on this plane 

> > 

such that pi G dToo{p{x),e) and p{x)p2, p{x)p^ be unit vectors lying on 

this plane such that P2;P3 £ dUoo{p{x)) and let \pi — P2\ < \pi — Psl- Then 

note that due to the chosen apertures of Too{p{x),e) and Uoo{p{x)) vectors 

> > I 

p{x)pi and p{x)p2 are perpendicular. It follows that H {x) is contained in 

U^{p{x)) and contains Ir{p{x)) (see Figure 2). 

Let li be the unit inner normal vector for H^. Let / be the point p{x) 

when the dimension d = 2 and let I be the (d — 2) dimensional hyper plane 

containing p(x) and perpendicular to u{p{x)) and n when d > 3. If d > 3 

and i'{p{x)) = let / be a fixed arbitrary {d — 2) dimensional hiperplane 

containing p{x) and perpendicular to v{p(x)). By rotating H-^ around I by 

some angle (3 we obtain the halfspace H'{x) such that x £ H'{x), 5jji(^^^{x) = 

^ir{p{x)){x) and Ir{p{x)) C H'{x) C U^{p{x)) (see Figure 2). Note that 

5ir{p(x)){x) = 8h'{x){x) < Sd{x) < \x-p{x)\ = 6h±(^,^){x) = 6u,(p{x)){x). 

It follows that by rotating H-^ around / by some angle, which is smaller or 
equal than P, we can obtain the halfspace H*[x) such that x G H*{x) and 

SH*{x)ix) = Sd{x). 

We also have 

Ir{pix))cH*{x)cU^ip{x)). 
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To simplify notation we will write 

fH{t,q)=rHit,{q,0,...,0),{q,0,...,0)), t,q>0, (2.21) 

where as before H = {(xi, . . . , Xd) G M'^ : xi > 0} is the half space. 

Our Proposition 2.9 below is a modification of Proposition 1.1 in [7]. We 
need this modification because we need to apply it to the function q 
fuO-iq) (see formula (2.21) for the definition of fn) which is not known to 
have the necessary estimates on its derivatives needed in the formulation of 
Proposition 1.1 in [7]. 

Let uj'^ = 7r''/2/r(s/2 + 1) and for < A; < d we recall that the the 
Minkowski content of a set E is defined by 

\{x e D : dist(2;,£') < r}\ 



M{E) = lim 



1, ^d-k 



0+ Wd-k^ 

We note that \i E <Z dD, the boundary of a Lipschitz domain, and E 
is closed, then M'^'^iE) = W^-^iE) [[13], Theorem 3.2.39], where W^'^ 
denotes the more familiar Hausdorff measure. Note also that ^u;^ = 1. 

Proposition 2.9. Let D C M*^ he a hounded Lipschitz domain. Suppose 
that f : (0, oo) —>■ R is continuous and satisfies f{r) < c(l A r^'^), r > 0, 
for some P > 1. Furthermore, suppose that for any < i?i < i?2 < / is 
Lipschitz on [i?i,i?2]- Then we have 

lim - / f (^^^] dx =n'^-\dD) r f{r)dr. (2.22) 
'7^0+ V Jd \ rj J Jo 

Proof. Let ipr]{r) = r]~^\{x G D : 6d{x) < ryr}|. We have (cf. proof of 
Proposition 1.1 in [ ]) 



?? ^ / f{^Dix)/rj)dx= / /(r)dV',,(r). 
Jd Jo 

So in order to prove the proposition we need to show that 

lim / f{r)diljJr) =rL'^'^{dD) fir) dr. (2.23) 
r?^o+ Jo Jo 

Let < i?i < ii2 < oo and r/ > be arbitrary. We claim that 

f{r)dil,^{r)<cRi, (2.24) 

JO 

oo 

/(r) #^(r) < cr]^-^ + cR\~'^ , (2.25) 
and that 

j-R2 rR2 
lim / f(r)d7pJr) =n'^-^(dD) fir) dr. (2.26) 

^^0+ Jr^ Jr^ 

Here and below constants c depend only on /, (3 and D and may change 
value from line to line. Note that (2. 24), (2.25), (2. 26) imply (2.23). 
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Proof of (2.24): Note that for any r] > the function r — > tpr}{f) is nonde- 
creasing and for any r/ > 0, r > we have ipr^ir) < cr. Hence 

/(r)#^,(r) < ||/||ooV'^(i?i) <ci?i| 



Proof of (2.25): Let M = sup^^jj 5£){x) (the inner radius of D). Since D 
is bounded, we certainly have M < oo. For i] > put = (M + l)/r]. 
Clearly for any r > we have tpni^) = il)rj{M^) = rj~^\D\. We have 



R2 



R2 



Mr, 



Integrating by parts this quantity equals 



cM-^Vr,(M^) - cR^'^ij^{R2) + c / r^/^- Vr,(r) dr 



R2 



which proves (2.25). 

Proof of (2.26): Set g{r) = \{x G D : 5d{x) < r}\. By definition we have 
lim^_o+ 9ir)/r = W^-^idD). Put 

9{s) 



r ^ dr 



R2 

,1-/3 



G{r) 



sup 

sG(0,r] 



Clearly G is nondecreasing and in fact G{r) | as r i 0+. Note that for 
r < i?2 we have 

'g{r]r) 



\ij^{r)-n''-\dD)r 



•qr 



W^-^idD) ) r 



< G{T]r)r < G{riR2)R2. 



Let P = {xo,xi, . . . , x„(p)} be a partition of the interval [Ri, i?2]- That is, 

Rl = Xq < Xi < . . . < Xn(p) = R2 

and let S{P,f,tprj), S{P,f) be the Riemann-Stieltjes sums 

n(P) 

S{P,f,'4'7^) = ^ f{xi){'ll;r,{Xi) - ll^rtiXi-l)), 
i=l 

n{P) 

S{PJ) = Y,f{xi){x,-x,.i). 
1=1 

Since / is Lipschitz on [Ri,R2] there exists a constant L > such that for 
ah ri,r2 G [Ri,R2] 

1/(^2) - f{n)\ < L\r2 - ri\. 
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This gives 



< 



\siP,f,^r,)-n''-\dD)s{p,f)\ 
f{xi){-i^^{xo) + n''-HdD)xo) 

n(P)-l 

+ E {f{x,)-f{x,+i)){i;r^{x,)-n''-\dD)3 

i=l 

+ f{XniP)){'4'vi^n{P)) - T-(-'^~\dD)Xr,(p)] 

n{P) 

2||/||ooG'(i?27?)i?2 + G{R2V)R2 - ^ 

i=l 



= G{R2v)R2mf\\oo + L{R2-Ri)). 
Hence for any partition P we have 

\S{P, f, - n''-HdD)SiP, f)\ < GiR2v)R2mf\\oo + LiR2 - Ri)). 
This impHes that 



R2 



f{r)di:^{r)-n'''\dD) I f{r)dr 



Ri 



R2 



< G{R2r])R2{2\\f\\oo+L{R2-Ri)), 

JRi 

which gives (2.26). □ 

3. Proof of main result 
Proof of Theorem 1.1. By (1.3) 

-t'^/° rD{t,x,x)dx = t'^/°' {pD{t,x,x)-p{t,x,x))dx 



D 



D 



t'^/'^ZD{t)-Gi\D\. 



Therefore to prove Theorem 1.1 we must show that for an arbitrary e > 
there exists a to > such that for any < t < toi 



f^'"" I rD{t,x,x)dx-C2n'^'\dD)t^''' 

JD 



< c{e)t 



l/a 



(3.1) 



where c(e) ^ as e — > 0. 

Fix < e < 1/4. We will use the (e,r)-good set G from Lemma 2.8 and 
corresponding set Q defined by (2.9). 

We need to estimate 

t'^/" / rD{t,x,x)dx. 



JD 

We split this integral into three sets 

Vi = {xe D\g : 6d{x) < s}, 
V2 = {xe Dng : 6d{x) < s} , 
V3 = {x£D: 5d{x) > s} , 



14 RODRIGO BANUELOS, TADEUSZ KULCZYCKI, AND BARTLOMIEJ SIUDEJA 

where s must be smaller than the sq given by Lemma 2.7. For small enough 
t we can take 



We can also assume that sq < r/4. 
By Lemmas 2.1, 2.7 and 2.8, 

t^/" / rDit,x,x)dx < c\Vi\ < C{dD)es = C{dD)./It^''^ , (3.2) 

where C{dD) is a constant depending on d, a and dD. 
Lemma 2.1 and the definition of the set 2?3 give 

t'^'"' / rD(t,x,x)dx < c / -T- Aldx 

< C / I - 1 a/e A 1 dx. 



Now we apply Proposition 2.9 with r/ = t^^" and /(r) = ^/er "^-"^ A 1. For 
small enough t this leads to 



t"^/" I rD{t, X, x)dx < cH'^-^{dD)t^/'^ / ^r-^-''^'^ A 1 dr 

(3-4) 



By (3.2) and (3.4), we have that 

t'^/" / rD{t,x,x)dx < c(e)t^/", 

where c(e) — > as e — > 0. 

It remains to consider the integral over I?2- Let x G I?2 C ^ and p{x) £ 
dD be such that x £ Tr{p{x),e). We also have 5d{x) < s. By (2.8), (2.9), 
(2.10), the fact that s < r/4 and e < 1/4 we have that \x—p{x)\ < 26d{x) < 
2s. Hence 5d{x) and \x — p{x)\ are comparable and x £ T2s{p{x),£). 

We now claim that on 1^2, r£){t,x,x) is comparable to rjj*(^^-^{t,x,x). Let 
Ir{p{x)), Ur{p{x)), H*{x) be defined by (2.17), (2.18), (2.20). Since 

Ir{v{x)) C H*{x) C U^Xpix)) 

and 

Ir{p{x)) CDC U^{p{x)), 

we have 

\rD{t,x,x) -rH*(x){t,x,x)\ < r/^(p(^))(t, x, x) - r[;c(p(^))(t, x, x). (3.5) 
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The next proposition asserts that for smah t, the difference ^/^(p) {t, x, x) — 
rjjci^p-^^t, X, x) is small. 

Proposition 3.1. For any p e M"^, any unit vector v{p) G M'^, e e (0, 1/4), 
t > 0, r > such that t^/"/i/e < r/A, s = t^/'^/y/e and x G V2s{p,£) we 
have 

0<r,,,(p)(t,x,x)-rt;c(p)(t,x,x) <c I [j^^^ 1 ^1 

The proof of this Proposition is fairly long and technical. In order not to 
interrupt the flow, we continue with the proof of Theorem 1.1 and return to 
the proof of the Proposition later. 

It follows from (3.5) that 

\rD{t,x,x) -rH*(cc){t,x,x)\dx < / r^^(p(^))(t, x, x) - rt7c(p(^.))(t, a;)dx 

(3.6) 

Let us first observe that 57,,(p(a;))(2;) < 5£){x) < 5[/^(p(2:))(a^) = \x — p{x)\. 
Notice also that by the definition of and Ir, and the fact that x £ 
T2s{pix),e), s < r/4 and e < 1/4, we have that x — p{x) and 57^(p(x))(3;) 
are comparable. That is, there exists a constant c G (0, 1) such that 
c\x — p{x)\ < Si^(p[x)){x)- Therefore c5d(x) < '5/^(p(2:))(x) and we conclude 
that c(^£)(x) and 57,,(p(a;))(3;) are comparable. This, together with (3.6) and 
Proposition 3.1 gives 

\rDit,x,x) - rH*(x){t,x,x)\dx 

2?2 




Once again we apply Proposition 2.9 with t] = t^l'^ and /(r) = r '^^^ A 1. 
For small enough t we get 



f^l"^ / \rD[t,x,x) -rii*(x){i,x,x)\dx 

r° 

< c(e^-"/2 V ^/I)n'^-\^D)t^/'' / 

Jo 

To complete the proof of (3.1), it remains to show that the quantity 



t"^/" / rH*(x){t,x,x)dx 

JV2 
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gives the second term C2'H'^"^{dD)t^/'^ in the asymptotics plus an error 
term of order c{e)t^^°^ . 
Recall that 

rH'{x){t,X,x) = fH{t,SH*{x){x)) = fH{t,SD{x)). 

Therefore 



= t""/" / fH{t, 6D{x))dx - t'^/" / fH{t, 5D{x))dx. 

Jd JV1UV3 
By the same arguments as in (3.2), (3.3), (3.4) we have that 

f"^" [ fHit,6Dix))dx<c{e)t'/'', 
JV1UV3 

where c(e) ^ as e ^ 0. 
Scaling now yields 

X ^Dix))dx = (1, ^) dx. 

By Lemmas 2.1 and 2.3 the function q fui^iO) satisfies the assumptions 
of Proposition 2.9. Hence we have that for small enough t, 

jjH{i,^-j^)dx-c,n'-\dD)t' 

This verifies (3.1) and finishes the proof of the Theorem 1.1. □ 

Proof of Proposition 3.1. We may assume that p = 0, z/(0) = (1, 0, ... , 0). 
To simplify notation let us define X = I^jii/a and U = 
Put 

I = {y -.yvi^) > e\y\}, 
U = {y:yiy{0) < -e\y\} 

T{0,e) = {y : y ■ iy{0) > Vl - e^\y\} 

By scaling (formula (2.3)) and Proposition 2.3 in [2] we have 

rir{t,x,x) - ruc{t,x,x) 

( X X \ (XX 



A/a 



\r-E-l^''^^) (tj < 1, Xirj) eU\I,pu(l- rj,X{Tj), 



X 



< ^E^/'^^'^"-^ (rj < 1, X(rj) ^U\I,p{l-Tj, X{ti), ^ 

Put u; = ^ G T^s/ti/c, {P, e) and 

R{w) = E'^in < 1, X{ti) eU\I,pil- TJ, X{ti), w)). 
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By the space-time Ikeda-Watanabe formula we have 

R{w) = / pi{q,w,y) v{y - z)p{l - q,z,w) dz dqdy. 
Ji Jo Ju\i 

Note that 

ZY\T = 5^(0,r/ti/°)U(C/"\/), 
Id. 

Therefore, 

Riw) < / Piiq,w,y) i^{y - z)p{l - q,z,w)dzdqdy 
J I Jo Jvxi 

+ [ [ Pi{Q,w,y) [ u{y - z)p{l - q,z,w)dzdqdy ^^''^^ 

JxJo J B'^{0,r/t^/°') 

= A{w) +B{w). 

In order to prove the proposition, it is sufficient to show that for any w E 

A{w) < cei-^/Vr"'"", (3-8) 
A{w) < ce^-'^l'^ (3.9) 

and that 

B{w) < c,/^\w\-'^~''+\ (3.10) 

Biw) < cVe. (3.11) 

Here we use the fact that 5/^(p(^))(x) ~ \x\ for x G T2s{p,£), s < r/4, e < 1/4. 
We first prove the easier inequahties involving B{w). Since s < r/4 and 

^ ^23/41/0 (p, e) we have 

SBc{o,r/t^/")iw) > cr/t^/". 
Therefore by (2.2) for q G (0, 1) and z G S=(0,r/ti/") we have 

/ , , \ d+a 

p{l -q,z,w) <c \ 



Substituting back into (3.7) leads to 

(l/a \ '^'^'^ 1 
/ / Pi{q,w,y) v{y-z)dzdqdy 
f I JjJo JB<=(0,r/ii/Q) 



' , , \ d+a 



<c(—] <c — <c{Vif\w\-''~''+^, 
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for any (3 < d + a. The inequalities (3.11) and (3.10) follow if we take /3 
equal d + a and 1 respectively. 

The proof of the inequalities (3.8) and (3.9) is quite technical. It will be 
divided into several steps. 

At first we will need the following auxiliary lemma. 

Lemma 3.2. For any e G (0, 1/4), w G r(0, e), M G (0, oo] we have 



The constant depends only on d,a,^. When M = oo, we understand 
B{0,M) = M.'^. 

In fact, we will use this lemma only in the cases when 7 = d, M = oo and 
when 7 = d — a, M = 1. We state it in this form two avoid repeating the 
proof twice. 

Proof. Let us introduce polar coordinates {p,ipi, . . . , tpd-i), with center at 
p = and principal axis z^(0) = (1,0, ...,0). There are some technical 
differences between the case d = 2 when ^pi G [0, 27r) and the case d > 2i 
when Lfi G [0,7r]. We will make calculations for the case d > 2>. The case 
d = 2 \s essentially the same, simply taking care of the restriction on the 
angle. 

Let G [0,7r] be the angle such that cos{'^e) = Note that 



and that 5i{z) = psm(ipi — c^^) for z £ I. 

LetVi = (C/^\/)nB(0,|w;|AM) andV2 = iU^\I)nB{0, M)nB%0,\w\AM) 
with the understanding that if \w\ > M, then V2 is empty. Note that 
|z — ~ \w\ for z €z Vi and — ~ \z\ for z G V2. We have 



The last inequality follows from the fact that for e G (0, 1/4) we have sin(7r — 
2ip^) ~ 2sin(7r/2 — ip^) = 2e, so vr — 2ip^ < ce. 




for 7 > d — a/2, 
/or < 7 < d - a/2. 




< 
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Similarly we have 

dz f dz 



Sv, 5^ \z - ^ fv2 5f^{z) \z 



J\w\AAI 

and the lemma follows. □ 
Now we will show that for w E r(0, e) we have 

/ Ki{w,z)dz <ce^-''/^, (3.12) 

where Kj(u!, z) is the Poisson kernel for /. For any z £ If^ \ I , let H{z) be 
a halfspace such that / C H(z), z £ {H{z)Y and 5i{z) = Sh(^z){z)- Recall 
that when H is a, half-space, yi £ H, y2 £ H'^, then (see e.g. (2.5) in [8]) 

KH{yi,y2) = ci , (3.13) 

% (y2)|yi - y2r 

where = T{d/2)T:~'^I'^-'^ sin(7ra/2). It follows that 



Kj{w,z)dz < / Kff(^z){w,z)dz 



< c\w 



a/2 



dz 



u-\i Sf^iz)\z -w[ 



This gives (3.12) by Lemma 3.2 for j = d and M = oo. 

We win now show (3.8). Note that by (2.2), for w £ r(0,e), z £ U" \ I 
and q £ [0, 1) we have 

1 — q c 
p[l — q, z, w) < C-. rr. — < -. — rr, — • 

\z — '^+" ~ kf K 

Using this, (2.7) and (3.7), we get 

M^) < I / / Pi{Q^w,y) I u{y - z)dzdqdy 



\yj\d+a 



Gi{w,y) I v{y- z)dzdy 



Now (3.8) follows from (3.12). 



Kj{w, z) dz. 
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Our next aim is to show (3.9). By (3.7) we have 

A{w) = [ [ Piiq,w,y) [ u{y - z)p{l- q,z,w)dzdqdy 
J I Jo Ju'=\i 

+ Piiq,w,y) i^iy - z)p{l - q,z,w)dzdqdy 

JlJl/2 J(C/<:\/)nB'=(0,l) 

+ / / PliQ,w,y) iy{y - z)p{l - q,z,w)dzdqdy 

JlJl/2 J {U'=\I)nB{0,l) 

= I + 11 + III. 

For q £ [0, 1/2], we have p{l — q,z,w) < c. Similarly for q £ [1/2, 1), w € 
r(0,e) and z £ (C/'=\/)n5^(0, 1), we have p{l- q, z,w) < c|t(; - z]"'^-" < c. 
Using this and (3.12) we obtain 

r roc p 

I + 11 < c / Pi{q,w,y) u{y-z)dzdqdy 
J I Jo Jvxi 

= c [ Ki{w,z)dz < ce^'"/2^ 

As for III, by [22, Theorem 1.6 and Corollary 1.7] we have pi{q,w,y) < 
cSf^iy) for g G [1/2, 1), y G I. Hence, 

m<c / I ' dy / pil-q,z,w)dqdz. (3.14) 

Thus for z£{U''\I)r\ B{0, 1) we have 

We also have 

1 roo ^ 

pi{l - q,z,w)dq < pi{q,z,w) 



'1/2 Jo \Z-W[ 

So by (3.14) we get 

dz 



III < c 



'{C/'=\/)nB(o,i) (^"^^(z) \z - w\'^-'^ 
Using Lemma 3.2 for 7 = d — a and M = 1, we finally arrive at (3.9). □ 
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